Looking for the Kondo cloud 
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We study the spatial distribution of the Kondo screening cloud. We find that in 3 D the spatial 
range of the cloud is limited to distances of the order of the Fermi wavelength due to angular 
dispersion effects. In lower dimensions this effect is less important and thus the Kondo length 
determines the localization of the cloud around the magnetic impurity. We clarify the role of the 
Kondo cloud in the impurity-impurity interactions: the hole coherence factor of the Kondo-doublct 
interaction is proportional to the amplitude of the Kondo cloud. We single out an experimental 
setup in which the Kondo cloud could be directly measured. 
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Introduction- When a metallic sample with magnetic im- 
purities is cooled below a certain energy an electronic 
correlated state develops, screening the impurity spin. 
That is the Kondo effect, probably one of the most well- 
studied phenomena in condensed matter physics both 
experimentally and theoretically. This effect is described 
by the Anderson Impurity Hamiltonian, whose properties 
in the different regimes have been studied by a variety of 
methods, ranging from Wilson's numerical renormaliza- 
tion group Q to "exact" Bethe ansatz results Q. Most of 
those results are covered in the book by Hewson, "The 
Kondo Problem to Heavy Fermions" [3| • 

Despite all this remarkable research effort there is a 
basic point that has eluded a clear theoretical description 
and experimental detection 0, It is the structure of 
the Kondo screening cloud and the role of the so-called 
Kondo length (£«•). 

The Kondo effect is characterized by a narrow res- 
onance of width 5k 1 pinned to the Fermi energy Ep. 
This resonance is intimately related to the formation 
of a many-body singlet state, comprised of the impu- 
rity states and a cloud formed from the metal band 
states that screens the impurity spin. Its spatial ex- 
tent is vital for the coupling between neighboring Kondo 
impurities 0]. Given the width and localization of the 
resonance it is straightforward to associate the length 
ii< = ^fEp/Sk ~ hvp/5k to it, where Xp (vp) is the 
Fermi wavelength (velocity). Although £k can reach al- 
most macroscopic values given that Ef/Sk — 1000, it 
has never been observed in experiments. 

The problem possed by such a large value for is 
rather disconcerting [H, El 0, S> 0- Even a very dilute 
system with few parts per millon of impurities has a typ- 
ical inter-impurity separation of about 100 lattice spac- 
ings, much smaller than Therefore each impurity has 
many other impurities inside its screening cloud and co- 
operative effects must arise. Nonetheless, the impurity 
resistivity, susceptibility, etc., are observed to be linear 
in impurity concentration up to very high concentrations, 



and those quantities fit the theoretical single impurity ex- 
pectations. 

In addition to the "classical" impurities in the 3 
D metal case, the Kondo effect has proven to be a 
very robust and general phenomenon, arising in 2 D, 
adatoms on metal surfaces[10j and quantum corral 11 
systems, and also in 1 D quantum dots-semiconductor 
heterostructures[12]|. In all these systems, and given the 
possibility to use Kondo quantum dots as quantum bits, 
it is thus relevant to clarify the structure of the Kondo 
cloud, its role in the interaction between neighboring im- 
purities, and the role of £k [H.[rj.[l3. 14|. 

In the following we solve this long standing puzzle. 

Theory. The Kondo effect physics is contained in the 
Anderson Impurity Hamiltonian. It is composed of a 
localized orbital (the impurity) of energy —Ed, measured 
from the Fermi level, with a strong internal Coulomb 
repulsion U, such that the energy cost to put a second 
electron in it (— Ed + U) is, in the Kondo regime, well 
above Ep. Thus double occupance of the impurity is 
effectively forbidden. This orbital interacts with a band 
of N extended states, characterized by a half-band width 
D, through an hybridization V, that allows an electron 
to be interchanged between the band and the localized 
orbital. In the Kondo regime U ^> Ed,D ^> 5 k the 
relevant Hamiltonian parameter is the effective Kondo 
coupling J = V 2 1 Ed and the Kondo energy is given by 
5k = D cxp (— 1/2 J n ), where J n = p J , p being the 
density of band states at the Fermi level. 

To analyze the Kondo cloud we use the Kondo singlet 
variational wave function (VWF) designed by Varma and 
YafetJUJ . This method was later used by Gunnarsson 
and Sch6nhammer[l(| to successfully explain the spectra 
of rare earth compounds, and has become known as the 
Variational l/Ns Expansion, Ref. [Hp. 223, where N$ is 
the degeneracy of the localized orbital, equal to 2 in the 
spin one half impurity case analyzed here. Actually, the 
expansion parameter can be shown to be J n /Ns- Typical 
values of J n are 0.1, for 5k — 0.01D, to 0.05, for which 
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8k — 0.0001-D. The results of this method are exact in 
the J n /Ns — > limit. The Kondo singlet VWF is given 
by 

i^> = i^>+E^ 6 L4i^>+E^ b t« c U F ) < C 1 ) 



for the singlet, obtaining 
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where \F) is the Fermi sea, i.e. the band filled up to 
the Fermi level. The second quantization operator &t 
creates a hole in the band, i.e. it removes an electron 
from the ka band-state below the Fermi level, c qa puts 
an electron in the band above the Fermi level and d\ puts 
an electron in the impurity. The k (q) sums are over hole 
(electron) excitations and a is the spin index. A figure 
showing these configurations can be seen in Ref.[4|. The 
variational amplitudes of the VWF are Z k = v / '(5 k + e&) , 
and Y kq = v Z k /(-E s + e k +e q ), e k (e q ) being the energy 
of the hole (electron) excitation in the band and v = 
V I SIN . The energy of the Kondo singlet is given by 



E$ = —Ed — St. 



Er 



(2) 



where the Kondo energy 5k comes from the resonance 
between the Z k configurations, which use the Fermi sea 
configuration as a nearly virtual bridge. Ej is the single- 
body impurity correction, equal to — J/2 for the level of 
approximation used in Eq.JT]). It comes from the inter- 
play between each Z ko configuration and their derived 
Yk oq configurations. Therefore Eq.([2]) ilustrates the ver- 
satility of this variational method: it gives both the non- 
perturbative 17[ Kondo term as well as the standard per- 
turbative ones. 

The square of the norm of \Sjc) is given by 
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and direct evaluation gives the relative total weight of 
the different set of configurations in the Kondo limit 
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where the usual (— i?s + e fc + e 9 ) ~ Ed approximation has 
been used in the denominator of Y kq . Thus practically 
all the weight of the VWF is on the Z k configurations, 
more precisely in the ones corresponding to holes near 
the Fermi level. One half of the weight comes from hole 
states between 5k of the Fermi level. To analyze the 
spatial distribution of the resonant hole we compute the 
band hole density operator, 



Ph( r ) = ^ E b L b k'a exp (i{k - k').r) , 

kk' (7 
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Eq.® can be rewritten as p K (r) = C ' K (6 K )\^> ' h {r)\ 2 
where 

= ^E e *" r ^ = ^ E ^ E elk ' r ■ (8) 

is the resonant hole wave function [3], and Ck com- 
prise the r-independent normalization factors. Given 
that Z k depends only on the energy of the hole exci- 
tation we separate the sum over k in its angular and 
modulus parts. The angular integral determines the prin- 
cipal dependence of on r and on the spatial dimen- 
sion D of the system. It gives cos kr, the order Bessel 
function J (fcr), and sin kr/kr, in 1, 2, and 3 dimen- 
sions respectively (i.e. ~ cos (kr — ira/2) / (kr) a where 
a=(D — l)/2). This behavior corresponds to that of 
radiation waves, and, although each k component is an 
extended state, it already provides a decrement of the 
hole cloud amplitude ~ /--(D- 1 )/ 2 . The dashed lines in 
FigO] correspond to those angular integrals evaluated at 
k = kp- 

In Fig[T] there are also plotted (full lines) the results of 
the whole sum, i.e. 'l'h(r). These sums can be analyt- 
ically done in 1 and 3 dimensions in terms of Cos/Sin- 
integral functions. The interference between the different 
terms of this last sum in |fc| generates true localization 
of the hole around the impurity, localization ruled by £k 
given that Z k ~ 1/ (5k + efc)- But, from a practical point 
of view, the amplitude of the Kondo cloud decays, in 2 
and 3 dimensions, in the order of a few times Xp due to 
its angular dispersion. For 5k — > 0, which corresponds 
to a loosely bound hole with k ~ kp , the ^(r) profile 
tends to that of the corresponding radiation wave, the 
dashed lines in Fig[T] 

The spatial extension of the Kondo cloud opens the 
possibility for the hole bounded to a given impurity to 
be captured by an other impurity located at a distance 
such that the orbital amplitude is still significant. Indeed, 
this is the that the hole coherence factor of the 

Kondo doublet interaction [3| is proportional to ^(r). 
This connection between the single impurity Kondo cloud 
and the recently discovered two-impurity Kondo-doublet 
interaction is not by chance. The Kondo-doublet interac- 
tion is mediated by the interchange of a hole between the 
impurities, i.e. the Kondo-doublet is formed by a hole 
which is partially screening two impurities. 

As stated in the introduction the Kondo effect has be- 
come relevant in many situations, each one with its own 
characteristics. 

In 3 D the angular dispersion effectively reduce the ex- 
tension of the Kondo cloud to few times Af, and £k does 



3 




FIG. 1: Kondo hole orbital profile for S K = 0.001D in 1, 
2, and 3 dimensions. The dashed lines correspond to the 
respective angular integrals cos(fcr), Jo(fcr) and sin(fcr)/(fcr) 
evaluated at k — kp- 

not play a significative role, in line with the observed 
experimental facts. Nevertheless, plots of r x h(r)\ 2 , 
with x = 2, 3, show that the hole density decays like 
1/r 2 up to very large values of r[20]. The charge den- 
sity around a 3 D Kondo impurity has been previously 
evaluated in Ref.0|,p.291 and [21[. Those calculations 
use mean field and/or perturbative solutions of the An- 
derson Impurity Hamiltonian and thus fails to capture 
the non-perturbative Kondo cloud. They found a 1/r 3 
dependence of the charge density that corresponds to the 
Friedel like oscillations, Sk independent, induced by the 
presence of a normal impurity in a metal. In Ref.[4| the 
higher order terms (~ 1/r 4 ) are attributed to the Kondo 
effect. 

In 2 D the angular dispersion effects are lower than in 
3 D and thus the Kondo cloud is effectively extended over 
a greater region. 2 D Kondo systems has been realized by 
the deposition of magnetic atoms over metallic surfaces 
with a strong density of surface states [loj. This situation 
is particulary interesting because electron density oscilla- 
tions can be straightforward measured by means of scan- 
ning tunnelling microscopy (STM), see Fig.l in Ref.[22| 
for Co adatoms on a Cu(lll) surface. In this setup there 



is a superposition of the Kondo cloud density oscillations 
and the Friedel screening charge oscillations [23]. But 
they can be discriminated by their r-dependence. The 
Freidel oscillations are given by 

A PF (r) ~ cos (2k F r + fo)/r u , (9) 

they not depend on Sk, and they normalize to minus the 
excess charge of the impurity, whereas that the Kondo 
cloud normalizes to one hole. 

In Fig[2]we plot r 2 \^ h {r)\ 2 for 2 D and different val- 
ues of Sk- The strong dependence of the Kondo cloud 
density with Sk can be seen. The charge Freidel oscil- 
lations (times r 2 ) are a uniform sinusoidal function in 2 
D. Instead, in Fig[2l for the Kondo cloud the decay rate 
is lower than 1/r 2 for a distance that scales with In 
fact, the maxima of r 2 |^ , ^(r)| 2 are at r/Xp ~ 5, 10, and 
20 for the cases plotted in FigJ2] 
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FIG. 2: 2 D Kondo cloud density profile (times r 2 ) for dif- 
ferent values of Sk- Its behavior is quite different from that 
of the charge Freidel oscillations (a uniform sinusoidal in this 
representation). The lower Sk, larger the density profile 
keeps decaying at a rate lower than 1/r 2 for longer distances. 

Another way to discriminate [3] between Friedel and 
Kondo contributions to the density oscillations is the fact 
that Kondo cloud disappears above the Kondo tempera- 
ture. Thus a comparative analysis above and below that 
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temperature will suffice to obtain the Kondo contribu- 
tion. 

In 1 D there are no angular dispersion effects, and thus 
the Kondo cloud has an appreciable amplitude over long 
distances, of the order of Given that in these systems 
formed by quantum dots/ wires in semiconductors het- 
erostructures A^ itself is very large, the door is open for 
the fabrication of very large integrated Kondo circuits. 
The Kondo-doublet interaction in them can be manipu- 
lated by means of tuneable spin-orbit couplings [HI and 
the spin-dot state measured by techniques like the one 
proposed in Ref. [25| . 

Higher order terms - The next configurations in the 
Kondo singlet VWF correspond, from the Yk q ones, to 
promote a new electron into the impurity, thus generat- 
ing a new hole in the band (i.e. b\, ^b^c^-dL, \F) ). 
These configurations are equivalent to a configura- 
tion plus a k'q hole-electron pair in the band. They con- 
tribute to the total energy of the singlet with a term of 
order p J 2 (= J n J). This is the order of the two impurity 
RKKY interaction[I3|. In fact the RKKY interaction is 
generate d by the annihilation of the eh pair at a nearby 
impurity 26| . The wave-function amplitude of this eh 
pair goes like l/r D , and its contribution to the density 
oscillations die off very quickly. The more excitations are 
present in a given term, more quickly it decays. 
Conclusions- We have found the Kondo cloud. Our anal- 
ysis of its spatial profile shows that in 3 D its range is ef- 
fectively reduced to a few times A^ due to the angular dis- 
tribution/dispersion of the orbital. The role of £k m de- 
termining the localization of the screening cloud around 
the Kondo impurity becomes more relevant in 2 and 1 D. 
We propose that the Kondo cloud can be directly mea- 
sured in a 2 D adatom setup, Kondo mirage effects have 
already been measured in a related geometry [27^. 

We also make clear the relation between the spatial 
extension of the Kondo cloud and the interaction with a 
nearby impurity. The amplitude of the cloud determines 
the possibility of the resonant hole to be captured by 
a neighbor impurity, and thus the hole coherence factor 
of the Kondo-doublet interaction is proportional to that 
amplitude. 

Beside the experimental possibilities opened by the 
knowledge of the spatial distribution of the Kondo cloud, 
further theoretical work is also in order. The Bethe 
ansatz solution starts from the proposal of an explicit 
wave function for the Kondo singlet, thus the details of 
the spatial distribution of the Kondo cloud could be also 
extracted from the ansatz. 
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